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NOMENCLATURE 


tip transient response-curve integrated areas, rad-sec 

tip forcing-frequency coefficient, nondim 

free-tip aerodynamic damping coefficient, N-m-sec/rad 

free-tip aerodynamic damping as a function of the nth harmonic 
component of the tip forcing frequency, N-m-sec/rad 

free-tip and controller-shaft mass centrifugal load, N 

tip lift coefficient curve slope, rad - ^ 

tip pitching-moment coefficient (nondim), positive noseup 

steady-state and nth harmonic component of the tip aerodynamic 
pitching moment (nondim), positive noseup 

unsteady nth-harmonic pitching-moment coefficient component that 
contributes to the tip aerodynamic damping, nondim 

unsteady nth-harmonic pitching-moment coefficient component that 
contributes to the tip aerodynamic spring, nondim 

tip inboard-edge chord length, m 

tip pitch-axis location, % of chord 

tip lift center offset from pitch axis, positive aft, m 

blade lead-lag hinge offset, m 

torsion-strap controller strap tensile load, N 

Theodorsen function (nth harmonic) 

second Theodorsen function (nth harmonic) 

free-tip and controller-shaft polar moment of inertia, kg-m 2 

tip-response reduced frequency, nondim 

constant-moment controller spring constant, N-m/rad 

free-tip aerodynamic spring constant, N-m/rad 

iii 
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free-tip aerodynamic spring as a function of the nth harmon 
component of the tip forcing frequency, N-m/rad 

K) equivalent-system spring constant, N-m/rad 

controller strap length, m 

free-tip and controller-shaft mass, kg 

controller moment output, N-m, positive noseup 

tip moment due to aerodynamic damping, N-m, positive noseup 

tip moment due to aerodynamic spring, N-m, positive noseup 

magnification factor, nondim 

air- jet plenum pressure, kPa 

controller strap radius, m 

tip forcing-frequency ratio, nondim 

rotor radius, m 

transient response envelope-area error estimate, nondim 
free-tip and controller output shaft radial c.g. location, m 
tip span, m 

total time interval of the free-tip transient response, sec 
controller coulomb damping moment, N-m 
tip downwash velocity, m/sec 
tip speed, m/sec 

weight of a free tip and controller shaft, N 
blade flap angle, deg 

tip transient response logarithmic decrement, nondim 

free-tip pitch angle, deg, positive noseup 

nth harmonic of the free-tip pitching response, deg 
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periodic forcing-function static displacement, deg 

torsion-strap controller strap pretwist, deg, positive noseup 

free-tip pitch angle with respect to inboard blade edge, deg, 
positive noseup 

controller strap slope, deg 

rotor rotational speed, rad/sec 

tip natural frequency, rad/sec, or uncertainty analysis accuracy 
estimate 

tip damped-response frequency, rad/sec 
tip forcing-frequency, rad/sec 

controller bearing surface coefficient of friction, nondim 

o 

atmospheric density, kg/m J 
blade azimuth, deg 

regression analysis, standard-error estimate 
tip transient response period, sec 
blade lead-lag angle, deg 
free-tip damping ratio, nondim 
free-tip/controller in-plane angle, deg 


v 



SUMMARY 


The aerodynamic response of passively oscillating tips appended to a model 
helicopter rotor was investigated during a whirl test. Tip responsiveness was found 
to meet free-tip rotor requirements. Experimental and analytical estimates of the 
free-tip aerodynamic spring, mechanical spring, and aerodynamic damping were calcu- 
lated and compared. The free tips were analytically demonstrated to be operating 
outside the tip resonant response region at full-scale tip speeds. Further, tip 
resonance was shown to be independent of tip speed, given the assumption that the 
tip forcing frequency is linearly dependent upon the rotor rotational speed. 


INTRODUCTION 


The free-tip rotor (FTR) is an innovative concept that has been proposed to 
meet two rotorcraft technology objectives: improve forward flight performance, and 

reduce rotor vibration loads. The FTR is similar to a conventional helicopter rotor 
for most of its span; however, the outboard 5-1051 is allowed to pitch with respect 
to the inboard blade. Figure 1 shows a model FTR. Tip motion is constrained by a 
mechanically applied moment that forces the tip to passively seek a constant lift 
level. Figure 2 shows the free-tip concept. The mechanical device that generates 
the applied moment that the tip aerodynamic and inertial moments oppose is called 
the controller. The FTR controller applies a nonzero moment with a negligible 
mechanical-spring rate, so as to minimize the variation in tip lift. 

Previous work published on the free-tip concept are references 1-10. Refer- 
ence 1 is the patent taken out on the FTR conceptual predecessor, the constant lift 
rotor. References 2 and 3 are the original FTR performance analyses. References 4 
and 5 contain the first FTR wind tunnel test data. Reference 6 is an early dynamic 
stability analysis of the FTR. References 7-10 are additional developmental studies 
of the FTR. 

The first wind tunnel test of an FTR model revealed that its controller did not 
perform satisfactorily (refs. 4 and 5); there was almost no tip response to aerody- 
namic excitation. Therefore, the FTR concept could not be properly evaluated during 
the test. A development program was initiated to design and demonstrate a success- 
ful controller mechanism (ref. 9). Several candidate controller configurations were 
subjected to rigid plane rotation and evaluated, and a controller design was 
selected for the FTR. Since the success of the FTR concept was contingent upon the 
responsiveness of the tip and controller mechanism combination to aerodynamic 
excitation, it was necessary to whirl-test the complete system to verify 
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responsiveness. The response characteristics of the system demonstrated the 
readiness of the FTR for wind tunnel testing. 

The purpose of this report is to present the results of the FTR whirl test. 
These results are in the form of tip transient-response measurements and empirical 
free-tip/controller spring and damping estimates. In addition to the experimental 
data presented, analytical expressions are derived for the free-tip aerodynamic 
spring, K ft , and aerodynamic damping, C A . Comparisons are made between the analyti- 
cal predictions and the experimental data. 


DESCRIPTION OF APPARATUS 


The model FTR used in this whirl test is the same model rotor used in the wind 
tunnel test described in references 4 and 5. The rotor modifications made for this 
test were a new controller and new tips that were increased from 5 to 1 0 % span. The 
model dimensions are given in table 1. Figure 3 depicts the FTR model in the test 
chamber . 

The free-tip controller design employed in this test is called the torsion- 
strap controller (TSC). The torque-output characteristics of the TSC mechanism are 
detailed in reference 9. The TSC derives its torque output from a twisted set of 
straps in tension, which is due to the centrifugal load of the tip mass, which is 
attached to the straps. The TSC design and the free-tip planform used in this test 
are shown in figure 4. 

For the rotor whirl test, tip aerodynamic excitation was provided by an air jet 
normal to the tip-path plane. Therefore a tip angle-of-attack change occurred once 
per revolution, which initiated a free-tip pitch-angle, 0, transient response. The 
air jet was located at an azimuth station of approximately 180°. The jet could be 
raised and lowered by a hydraulic jack, and had an operating plenum-pressure range 
of 0 < AP < 41 kPa. The test-chamber corners were located at i|> = 45°, 135°, 225°, 
and 315°. The interference effects of the test-chamber corners were found to have 
significant effects on the free-tip responses. 

A Hall-effect transducer measured the free tips' oscillations in response to 
the air jet, and was used to measure the pitch displacement, A9, of the free tips 
with respect to the inboard blades. The transducer measured the magnetic flux of a 
magnet mounted on a free tip; as magnetic flux density is directly related to the 
distance from the magnet, the output voltage could therefore be calibrated for free- 
tip pitch-angle displacements. The response curves were sampled and averaged by a 
computer and reduced to Fourier coefficients. 
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FREE-TIP WHIRL-TEST RESULTS 


The FTR concept requires the free tips to be highly responsive to aerodynamic 
forcing (u> >> 1/rev) and to have a low controller spring constant relative to the 
aerodynamic spring (K << K^). The purpose of the transient response data acquired 
during the whirl test was to verify that the free-tip/controller frequency and 
spring requirements were met. Additionally, analytical predictions of free-tip 8 
response have not been verified because of a lack of experimental data. 

Before proceeding to the development of the experimental data-reduction algo- 
rithms and the free-tip aerodynamic theory necessary for data correlation and the 
verification of the tip/controller design response characteristics, it is appro- 
priate to make a number of qualitative observations concerning free-tip 0 time 
histories and FTR performance data acquired during the whirl test. 


Free-Tip Time Histories 

Free-tip 0 time histories were measured by a Hall-effect transducer. Each 
time-history curve is the average of 20 tip time-history cycles. The individual 
0 time-history cycles were not retained by the data-acquisition system. All time 
histories were reduced by Fourier analysis to a set of experimental harmonic coeffi- 
cients (summarized in table 2). A typical time history that includes the free-tip 
response to air-jet excitation is shown in figure 5. 

Figures 6(a)-(c) are typical air-jet velocity distributions for various plenum 
pressures and displacements from the tip-path plane. The magnitude of the free-tip 
impulsive deflections from the steady-state pitch angle was dependent upon these two 
air- jet parameters. The air- jet velocity profile is assumed in the transient 
response data-reduction analysis to approximate a step-function excitation of the 
free tips. 

Figure 7 is a transient response curve that has been digitally filtered to 
remove the steady state, the 1/rev, and on a conditional basis, the 2/rev harmonics 
of the tip pitch angle. The steady-state 0 offset was removed to analyze the 
transient response. The 1/rev was eliminated to remove lateral and longitudinal 
cyclic influences on tip behavior. The 2/rev harmonic was sometimes filtered out 
according to a minimization of error strategy for the and coefficient 
regression. The 2/rev filtering was performed in an attempt to reduce the test 
chamber's corner-interfShence effects on the tip response data. Because a set of 
open windows were located along one wall of the test chamber, there was probably 
minimal aerodynamic interference from two of the test chamber's corners. For this 
reason the 2/rev was filtered out of the transient response data instead of the 
4/rev. The initial deflection of the figure 7 transient response curve is located 
at a rotor azimuth position (i|> « 190°) corresponding to the free tip exiting the air 
jet, where 0 « 0. 
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The transient response of a free tip subjected to a step-function excitation 
should theoretically be the decaying sinusoid of a single degree of freedom, linear, 
combined coulomb and viscous underdamped system. A key reason for the nonlinear 
disagreement between the theoretical response and the experimental data may be the 
aerodynamic influence of the test chamber corners — in particular, the first corner 
after the air jet probably initiates a tip Aa change. The corner creates a recir- 
culation zone that provides a localized upwash that interferes with tip response. 

It was therefore necessary, in developing data-reduction algorithms, to accommodate 
this response-curve distortion in the aerodynamic spring and damping calculations. 


Whirl-Test Rotor Performance 

During the FTR whirl test, rotor-performance data was acquired at very low 
C T /o. Figure 8 is a f igure-of-merit curve, and figure 9 is a power-versus-thrust 
curve for the model FTR tested. Rigid-tip baseline information was not acquired. 
Low-thrust performance curves were obtained for two different controller pretwists 
(output moments) applied to the free tips. Rotor performance was not significantly 
affected by the free-tip control moment for the limited data set acquired. 


FREE-TIP RESPONSE EXPERIMENTAL DATA-REDUCTION 


To gain a complete understanding of the free-tip/controller response character- 
istics, it is necessary to determine values for the coefficients K, T f , K^, C A , and 
I, which are, respectively, the controller mechanical spring and damping, the tip 
aerodynamic spring and damping, and the tip/controller polar moment of inertia. In 
the analytical results section, expressions will be derived for these response 
coefficients. In this section, data-reduction methods are developed that will later 
be used to calculate K, K A , and C A from experimental data for comparison with 
analytical predictions. Unfortunately, a satisfactory method for empirically deter- 
mining Tf was not developed. The remaining coefficient, I, was estimated by 
simple analysis of the tip density distribution. 

The mechanical spring, K, will be obtained from regression analysis of control- 
ler pull-test data. Three different methods will then be developed that attempt to 
estimate the aerodynamic spring and damping. All three data-reduction methods use 
the free-tip, air-jet response data to calculate K A and C A . 


Controller Experimental Coefficients 

Empirical estimates of the controller coulomb damping are unavailable. The 
data-reduction methods used could not satisfactorily isolate the Tf contribution 
from the total damping of the free-tip response curve. 
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Calibrations were performed on the controller that allowed an estimate of the 
spring constant to be made. Two check-load calibrations were made. First, a pull- 
test rig was assembled such that a tensile load (F = 7*100 N, 100£ rpm) could be 
applied to the TSC, while letting the strap angle be indexed through its allowable 
range. Second, the TSC controller was set to a specified angle (e = 70°) and the 
strap tensile load was varied within its structural limits. The calibration curves 
are shown in figures 10(a) and 10(b). 

The controller spring constant, K, can be estimated by these two curves by 
least-squares regression (ref. 9). The controller control moment can be represented 
by a semiempirical expression where the constant, A, is determined by regression 
analysis : 


M = A0F 
c 


Therefore 


K = AF = Amfl 2 R 

m 

An expression for A can be derived by noting the above M Q expression and equa- 
tion (11) in the analytical-results section. 


A = 


r_ 


When the two controller calibrations were performed, it was noted that there 
was a 1 0% difference between the experimental estimate for A_ and the calculated 
analytical value 
cally derived A 
from the empirical value for 
(800 rpm) . 


Experimentally, A was found to be 3.5x10 -4 m, and the analyti- 
is equal to 3.9x10"^ m. Estimating the controller spring constant 
A results in K = 2.59 N-m at V T = 213 m/sec 


A typical analytically derived value of the aerodynamic spring is 
Kft n = 53 N-m, as will be demonstrated in the analytical- and experimental-response 
comparison section. The subscript n denotes the number of the K^ harmonics 
used, which was 6. This choice of harmonic number will be explained in the experi- 
mental- and analytical-results correlation section. Therefore, with CF aligned 
with the controller straps, K contributes only approximately 5 % of Kg. 

It would be extremely difficult to isolate the actual K contribution to the 
K E derived from the transient response data, because of the difficulty in defining 
the in-plane and out-of-plane free-tip/controller angular orientation. For although 
pull-test calibration data were obtained, the strap alignment with respect to the 
tip centrifugal force is not known and, as will be shown in the analytical-results 
section, a value for the controller spring constant cannot be estimated from the 
whirl-test data. Therefore, since the K values are relatively small, the experi- 
mental Kg estimates will be compared directly to analytical predictions of K^. 
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Aerodynamic-Response-Coefficient Algorithms 

Three data-reduction algorithms were developed to calculate the free-tip 
response coefficients from the experimental data. These algorithms were a 
logarithmic-decrement-averaging method, a differential-equation, least-squares- 
regression method, and a linear-decay-averaging method. Table 3 shows each method's 
coefficient-analysis capability. The differential-equation, least-squares regres- 
sion approach is theoretically the most versatile of all the methods in that it 
should be capable of estimating all of the free-tip response coefficients. In 
comparison, the other two data-reduction methods were used to analyze the experi- 
mental data in terms of simplified models of the tip dynamic-pitching equation, 
resulting in a partial coefficient set. For example, the logarithmic decrement 
method assumes the tip response is only viscous damped, whereas the linear decay 
method assumes only coulomb damping. 

The relative contributions of coulomb and viscous damping to free-tip damping 
is unknown because T*> cannot be experimentally determined. Further, using an 
analytical T^ prediction to determine the contribution of coulomb damping is 
impossible, as will be demonstrated in the analytical-results section. However, the 
accuracy of each of the above simplified models must be gaged by the accuracy of the 
curves generated from the coefficient estimates. For example, a poor comparison of 
an experimental 0 time history with a response curve derived from a linear decay 
model would indicate that the tip/controller response is not predominately coulomb- 
damped. 

The coefficient-regression results indicate that the logarithmic-decrement 
method is superior to the other two methods. Therefore, the free-tip/controller 
damping is predominately viscous. The logarithmic-decrement method will be dis- 
cussed below. The remaining two methods will be discussed in Appendices A and B. 

The logarithmic-decrement method is used to determine Kg and C fl . T^ cannot 
be estimated by this analysis. The method defines a "local" set of dynamic coeffi- 
cients, based upon each cycle's amplitude decay, and estimates the tip-response 
coefficients by averaging all the local estimates. Equations (1) through (7) are 
the applicable equations for this method. 

The transient response amplitude decay is represented by the parameter 6 



where 0j and 0j + ^ are the adjacent, transient, angular-displacement peaks for 
the jth and (j+1)th cycles. If the experimental response curve demonstrated no 
coulomb damping, nonlinear behavior, or data scatter, the parameter need be evalu- 
ated only once to estimate the system dynamic characteristics. However, to minimize 
the influence of these effects, the experimental data are evaluated for each cycle 
to obtain several 6 estimates. 
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Next, given equation (1), reference 11 defines the damping ratio as 
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( 2ir )‘ 
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The remaining dynamic parameters can be determined from the following 
equations. 


( 2 ) 


to 


2 IT 

d . - t . 
1 1 


(3) 


where to d is the damped frequency, and is each cycle's damped-response period. 


to 


to . = 
l 


d. 

l 


V 


(4) 


1 - c 


d. 

i 


The natural frequency, to, is related to the damped frequency, to^, by means of the 
damping ratio. 

Therefore, K E and are estimated locally to be 


C A. = 2c d. Iu i 


(5) 


and 


K e = IOK 
1 


( 6 ) 


where I is the free tip and the controller polar moment of inertia. 


Finally, the effective response coefficients are assumed to be averages of the 
local, per cycle, coefficients. 

N 


c r » Z c a. 

i= 1 1 

N 


(7) 






i=1 
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Regression-Analysis Error Estimates 


As was noted in the whirl-test results section, the free-tip time histories 
were digitally filtered prior to data-reduction . The steady-state and 1/rev har- 
monic terms were excluded from the Fourier series representing the tip transient 
responses. The 2/rev harmonic terms were sometimes also removed from the tip- 
response Fourier series. The 2/rev removal was performed in an attempt to reduce 
nonlinear distortions to the tip transient responses. These response distortions 
are not representative of the theoretical linear-response curves and are perhaps due 
to the interference effects from the test-chamber corners, of which were not suc- 
cessfully filtered completely out. Sets of coefficients were obtained from both 
versions of the transient response curve: one curve retaining the 2/rev, and the 

other with it removed. The results of both coefficient regression analyses were 
compared to the filtered experimental data and error estimates were calculated. The 
Fourier series, either with or without the 2/rev, which produced the smallest error 
estimate was assumed to have the more accurate coefficients. Therefore, the accu- 
racy of the free-tip/controller spring and damping coefficients determined how much 
digital filtering was applied to the data. 

Arriving at suitable error expressions was, therefore, an important data- 
reduction concern. The following section discusses the two error estimates that 
were used on the whirl-test free-tip data. 

The first error estimate used on the whirl-test data was the standard error, 
which gave a point-by-point comparison of the experimental and curve-fit tip tran- 
sient response. The second error estimate is concerned with free-tip gross-response 
behavior. It is perhaps more suitable as an exactness-of-f it measure than the 
standard error estimate, considering the test-chamber interference effects are not 
modeled, and therefore tend to invalidate a point-by-point comparison. This error 
estimate will be referred to as the envelope-area error estimate. 

Using the envelope-area error measurement, a comparison is made between the 
experimental and curve-fitted response-curve envelopes. The assumption behind this 
error estimate is that equivalent envelope areas should indicate approximately 
equivalent envelope coordinates. Refer to figures 11(a)— (d). Figures 11(a) 
and 11(b) define the response-curve envelope and the envelope area respectively. 
Figures 11(c) and 11(d) demonstrate in principle an envelope-area comparison between 
the curve-fit response and the experimental response. While interference effects 
will negate a point-by-point correlation, their impact on parameter comparisons such 
as the envelope area should be less severe. The applicable equations for the 
envelope-area error estimate are shown below. 

The experimental envelope area is approximated by a first order, trapezoidal- 
rule numerical integration. 


exp 


■I 

i=1 


ASS 2 (e i+1 


8 i )(fc i + 1 


V 


( 8 ) 
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where 0^ and are adjacent experimental transient-peak amplitudes for each 

cycle, and t^ and are the corresponding response times. 

The curve-fit envelope area can be analytically estimated by integration of the 
exponential decay term in the underdamped-vibration general solution. This results 
in the expression 


where 



T = e 


(-3irc d )/(2 


V 




(9) 


Finally, a nondimens ional parameter was derived that was used to compare not only 
the curve-fit results for a single test point, but also the results for other com- 
patible test points. This parameter, R £ , eliminated the dependence of the error 
estimate on the free-tip initial deflection, which varied from one test point to 
another. Further, this particular parameter was used in weighted averages of the 
curve-fit coefficients. Coefficient weighting is appropriate in obtaining coeffi- 
cient averages in that the relative accuracies of the regression analysis will 
differ significantly between data sets. The weighted averaging will be discussed in 
the analytical and experimental coefficient-comparison section. R £ is defined as 


where 


ABS( A 


R„ = 


exp 


A calc^ 


norm 


( 10 ) 


A = e t 
norm o 

is the normalizing scale, and t in equations (9) and (10) is the total time inter- 
val of the tip transient response. 


FREE-TIP-RESPONSE ANALYTICAL RESULTS 


A free tip is modeled as a single degree of freedom system with its controller 
and aerodynamic springs, and dampers, in parallel. In general, the free-tip 
response equation has the form 

10 + C ft 0 + (K fl - K)0 + sgn(0)T f = M(t) 
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where I is dependent upon the chordwise mass distribution of the tip and control- 
ler, and are dependent on the tip pitching-moment coefficient, K and are 
functions of the controller kinematics and applied loads, and M(t) is the time- 
dependent tip aerodynamic and centrifugal forcing moment. The function sgn(e) 
returns the sign of 8. 

The differential equation can be simplified by noting that there is no continu- 
ous external unsteady forcing being applied to the free tips during the whirl test 
(i.e., tip speed and the flap and lead-lag angle are all steady). The air jet does 
impose a 1/rev, short duration forcing of the tips, but for most of the rotor azi- 
muth range, the tips are free of unsteady forcing. Therefore, the term M(t) 
becomes negligible when the tips have exited the air jet, and the response equation 
becomes 


1 9 + C^0 + (K^ - K) 0 + sgn( 0 )Tj. = 0 

The objectives of the following analytical derivations are 1) obtain expres- 
sions for the dynamic parameters included in the whirl-test, free-tip response 
equation (i.e., C^, K^, K, and T^ ) ; 2) demonstrate the validity of a simplification 
in the response equation that allows comparison between analytical predictions and 
experimental data (i.e., K << and Kg a K ft ) ; 3) demonstrate analytically that 
the free-tip pitching motion is stable and nonresonant. 


Derivation of Controller Response-Coefficients 

The following linearized expression for the torque output of the TSC is pre- 
sented in reference 9. The effect of flapping and lead-lag are not incorporated 
into the model. 


M 


c 



- A0 ) 


(ID 


where F is equal to CF. 

A TSC control-moment expression incorporating out-of-plane motion (flapping and 
lead-lag) will now be derived. The TSC derives its torque from a component of the 
strap tensile load, which is a component of the CF of the controller and tip 
mass. Since the controller out-of-plane motion will reduce the component of the 
CF acting through the TSC strap spanwise axis the torque output will also be 
reduced. Therefore, the condition for the maximum control moment occurs when 
8 = 0 and c = 0. 

Refer to figure 12 for the freebody diagram of the effect of 8 and t on the 
tensile load. It can be seen that the strap tensile load, F, is expressed by 

F = CF cos 8 cos £* 
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where 


= sin 


’IV 

R 

L m 


sin( 180° - c) 
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Thus, equation (11) becomes 


M * 
c 




(e pT - A0) 


which has the form 


M c - K(e pT - A0) 

and so the controller spring is therefore 


K = CF cos e cos c* yfy (12a) 

A semi-empirical approach using equation (12a) and a check-load calibration is 
found to be necessary to accurately predict the TSC spring constant. This will be 
discussed in the analytical and experimental coefficient comparison section. The 
semiempirical estimate for K will confirm the relatively insignificant 
{* 5 % of K ft ) contribution of K to Kg. 

Another controller parameter that contributes to the dynamic (and static) 
behavior of a FTR free tip is the T p term, the mechanism's friction-moment compo- 
nent. Reference 9 provides an expression for T f for the special case of rigid- 
plane rotation: 



yr V(w cos 6 + CF sin. 8) + (CF sin c#)‘ 


This expression incorporates 6 and c contributions to T p . However, as a 
rigid rotating plane is assumed, the inertial loading on the controller bearings 
caused by the derivatives d 2 B/dt 2 and d 2 c#/dt 2 is not considered. An extension of 
the above T p model will now be made to include the out-of-plane accelerations. 


Figure 13 is a freebody diagram of the TSC showing its bearing contact forces, 
including inertial loads. In defining the inertial forces it is assumed that the 
tip and controller output shaft CG is located on the shaft axis. The inertial 
forces are additive to the steady-state bearing contact forces. Therefore, the 
modified friction equation becomes 


= yr 


)[ v 


cos 8 + CF sin 8 + mR 


m 


9 • (■ 


CF sin c* + mR 


d 2 c»\ 
m dt 2 ) 


(12b) 


11 


The controller out-of-plane angular orientation and acceleration components are 
very difficult to determine. Additionally, y is not easily determined. Therefore, 
the controller T f could not be analytically predicted for the whirl-test case. A 
measurement of the controller friction coefficient and controller accelerations are 
required before an accurate T^ estimate can be made. 


Derivation of Aerodynamic Damping Coefficient 

In forward flight, the free tip is subjected to a periodic, external forcing 
function consisting of tip relative-velocity changes, varying local inflow, and 
blade flapping. A rotor in hover, or being whirled, is generally not subjected to a 
tip forcing-function. For the free-tip whirl-test, the tips are subjected to forc- 
ing once per revolution by sudden angle-of-attack changes caused by an air-jet 
blowing through the tip-path plane. Considered on a per cycle basis, the air-jet 
excitation approximates a step function input to the free-tip spring and damper 
system. Given the relative simplicity of the whirl-test free-tip forcing, it is 
possible to analytically derive C fl and K^. 

An analytical expression will now be derived for the aerodynamic damping coef- 
ficient for a whirl-tested free-tip. There are two important assumptions for the 
following aerodynamic damping and spring analysis: steady inflow and negligible 

periodic out-of-plane blade motion. For a rotor in hover, or being whirl-tested, 
these assumptions are valid. 

Equation (13), derived in reference 3, describes the free-tip aerodynamic 
moment 


M„ = \ C c 2 ARp(ftR) 2 
A 2 m o 

with the aerodynamic moment consisting of two components 


(13) 



+ M d 


where Mg is the moment due to aerodynamic spring, and is the aerodynamic 

damping moment. 

The contribution to caused by damping is 

M d = I C md c o 4R » (aR » 2 (lit) 

where C md incorporates all the terms of C m that are constant multipliers of 
de/dt. 


Reference 10 demonstrates that C m can be represented by a linear combination 
of harmonic components for unsteady flow about the free tip. Therefore, the damping 
moment contribution for each harmonic can be written as 


12 


( 15 ) 




ARvfc . 
t md 


n 



where V fc = nR and the subscript n denotes the harmonic number. 


Reference 10 derived equations (I6(a)-(d)) for the harmonic components of C m , 
for the two-dimensional (2-D) case. These equations are based on Theodorsen 
unsteady potential-flow analysis. It is necessary to apply corrections to the 2-D 
inviscid and K fl predictions or they will significantly differ from the experi- 
mental spring and damping values. The possible corrections might be in the form of 
modifications to the 2-D, unsteady pitching -moment coefficient equations or substi- 
tution of new Cfc values obtained from computational fluid dynamic solutions or 
from experimental data. Data from a low-speed, wind-tunnel semi-span test with 
indexed tips was used to calculate Cj^ values. The geometry of one of the tips 
studied during this test was the same as the whirl-test free tips. Therefore, an 
attempt to correct the damping predictions for the free tips was made by using these 
empirical tip semi-span values for in the 2-D, potential-flow C m expression. 

The steady-state pitching-moment coefficient is obtained from 


: = <=. i(° - * f v ♦ o (i - o ) ^1 , °o / _ i) 

m o a a\\ pa 4 /L° V t °\ 4 P a / V tJ 4 \P a V v t J 


(16a) 


The harmonic unsteady-moment coefficient is obtained from 


c = c„ 

m n 1 


o ((°pa - *){£» - (? - G pa) 

[. (I - °pa) -„)} * r {( G pa - i) (£ * <“„) 

t t ' 


where the reduced frequency, k, is 


k = 


c nil c n 
o o 

2V^ = 2R 


( 1 6b ) 


( 16c) 


and the functions, F n and G n , are defined by the complex Bessel-function expression 
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( l6d) 


F + iG = 


J^k) - iY 1 (k) 


J 1 (k) + Y o (k) - i [ Y 1 ( k ) - J Q (k)] 


Damping coefficients are defined as constant multipliers to d0 n /dt. There- 
fore, the 0 time-derivative coefficients are separated from the above equations 
(I6a-d) and are incorporated into the C md expression. And so 


C . = C„ 

md l 
n 


Equation (15) can be rearranged as 


.{('pa - i)[-o(l - 4] ^ (°pa - }) 


(17) 


M 


d0 /dt 
n 


where 


C A = 5 “‘o aBV ? C md 

n n 


( 18 ) 


Substituting equation (17) into (18), and simplifying, results in the following 
relation 


C = 1 o arc 

n 


o [ C 0 (l - °pa)(l - c pa) F „ V T * S °o (f - c pa) °n 

M (J - «J»t] 


(19) 


All that remains to estimate the free-tip aerodynamic damping are the explicit 
expressions for both F n and G n . By manipulating equation (I6d) and separating the 
real and imaginary terms into two relationships, the two functions are found to be 

J 1 (k)[J 1 (k) + Y q ( k) ] + Y^kHY^k) - J Q (k)] 


[J^k) + Y o (k)] 2 + (Y 1 (k) - J Q (k)] 2 


J 1 (k) [Y 1 (k) - J Q (k) ] - Y 1 (k)[J 1 (k) + Y Q (k)] 
[7, (k) - J Q (k)] 2 + [J^k) + Y Q (k)] 2 


( 20 ) 


where J Q (k), J-|(k), Y Q (k), and Y^(k) are Bessel functions of the first and second 
kind of order zero and one, dependent upon the reduced frequency, k. 

Equations (16c), (19), and (20) describe a complete solution for CA n > 
free-tip aerodynamic damping. An uncertainty analysis of the aerodynamic-damping 
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equations is included in Appendix C. The predictions discussed in the 

analytical- and experimental-results comparison section will be accompanied by 
damping uncertainty estimates based upon the analysis in the appendix. 


Derivation of Aerodynamic Spring Constant 


Similarly, an equation for the aerodynamic spring can be estimated from coeffi- 
cients taken from equations (l6a-d), and substituted into equation (13). 


M 

s 



AR^C 
T ms 


n 


0 

n 


where C mSn incorporates terms from equations (16 a-d) that are constant multi- 
pliers of e n . 

The spring is then 

M 

K A - 0 
n n 

Or, therefore 

K a =ipc2 a R V 2c ms (21) 

n n 

Where from equation ( 1 6b ) 

C ms n ” C * a ( C pa " ^)j^n " C o(4 " °pa) R ] (22) 

Substituting equation (22) into (21) results in 

\ - I 0 (if - C pa)- C H a [ F n - °o (l - c pa) ¥] <23) 

An uncertainty analysis of the aerodynamic-spring-constant equations is dis- 
cussed in Appendix B. 

The coefficients for the free-tip response equation have all been derived. An 
analysis defining the conditions for free-tip resonance is next described. 


Free-Tip Pitching Resonance 

An important objective of FTR research is avoiding resonant tip pitching oscil- 
lations. Given the derived approximate expressions for and K^, two important 
results can be derived for a free-tip subjected to periodic forcing. 
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The first result is an equation relating tip geometry and aerodynamics to the 
tip oscillation amplitude, expressed in terms of the magnification factor. The 
magnification factor (MF) is the ratio of a tip's actual oscillation amplitude to 
the tip displacement that would result from a statically applied moment equal in 
magnitude to the maximum periodic forcing. 

The second result is the mathematical demonstration that tip resonance is 
dependent upon tip geometry and inertia but independent of tip speed, given the 
assumption that the aerodynamic forcing frequency is a linear function of ft, which 
is justifiable because helicopter rotor aerodynamics and dynamics are often n/rev 
dependent. 


Reference 1 1 defines MF to be 


where 


MF 


0 _ _ 1 

0 f 272 ~ 72 

y(1 - r f ) + (2e d r f ) 


f 

s — (forcing frequency ratio) 
C A 

t d = 2^7 ( dam P in g ratio) 


0) 



(tip natural frequency) 


(24) 


The free-tip pitching-resonance analysis presented in this paper will be limited to 
pure sinusoidal forcing. The more complex case of general periodic forcing, like 
that of the whirl-test air jet or forward-flight rotor aerodynamics, will not be 
considered here. Additionally, it should be noted that tip-resonance analysis for a 
rotor in forward flight would require new aerodynamic spring and damping equations 
that include the effect of blade motion and unsteady inflow. However, such an 
extension could be made by representing the general periodic forcing by a Fourier 
series. The analysis will be further simplified by ignoring controller spring and 
damping contributions to the tip response. 

The Vp and c d terms in equation (24) need to be derived for an FTR free 
tip. The forcing-frequency ratio, r^, will be derived first. 

Let the forcing frequency be a linear function of ft, or 


Wj, = bn 


In terms of tip speed 
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0) 


f = 



(Note that if b is an integer, then the tip forcing frequency is n/rev.) 


To complete the definition of rp it is necessary to derive an expression for 
the natural frequency of the free tip. To do so, the expressions for and 
must be substituted into the natural frequency definition. Note that the dependence 
of the aerodynamic spring on tip speed is parabolic while the damping coefficient is 
linear. 


where 



(25a) 



B and D are constants that are obtained from equations (19) and (23). Note that 
G n and F n are terms that are comprised of Bessel functions that are dependent only 
on k, the reduced frequency, and not on tip speed. The rotor rotational-speed 
terms are canceled out in the k expression (see equation (16c)). 

Therefore, in terms of the above defined constants 


<D 




V 


T 


and so, from the previous work 



Therefore, r^ is independent of tip speed. 


( 26 ) 


An expression will now be derived for c d , the damping ratio. Substituting 
the and K A functional expressions into the definition results in the 

relation 
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( 27 ) 


C d " 2Iu> “ 
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_D 

^ 2/lB 

The damping ratio is also independent of tip speed. 

The MF has been shown by equations ( 24 ) — ( 27 ) to be independent of rpm. There- 
fore, as long as the assumptions < 0 - » bn and K << hold true, then free-tip 
resonance (MF > 1 and r f = 1) will not be initiated by rpm changes. Additionally, 
it should be noted that the solution of equation (19) cannot result in a negative 
estimate for C^. Therefore, the damping analysis indicates that free-tip oscilla- 
tions will be stable in hover. 


COMPARISONS BETWEEN ANALYTICAL AND EXPERIMENTAL RESPONSE COEFFICIENTS 



Application of the coefficient regression algorithms to the experimental data 
revealed that both the coulomb only, the combined coulomb and viscous damping analy- 
sis options of the least squares method, and the linear-decay averaging method 
resulted in erroneous coefficient values. The results obtained by these methods 
were physically unrealistic, i.e., application of either method would frequently 
result in one of the coefficients being negative. While free tip responses were 
obviously underdamped and not self-excited, these results were considered invalid 
and were therefore not included in this paper. 

It was found necessary for almost all curve-fits to exclude a portion of each 
response curve so as to avoid test-chamber corner aerodynamic-interference 
effects. The azimuthal location of the air jet was approximately tp = 190°. There- 
fore, the curve fitting was often limited to data having azimuths greater than 
4) = 240°. To better demonstrate the validity of the curve fits, the resulting 
curves were extrapolated into the excluded azimuth range. How well these two sets 
of curves match is an indication of both the magnitude of the corner-interference 
effects and the quality of the curve-fit. A good curve-fit based on experimental 
response-curve data with minimal interference effects has a smooth transition 
between the two curve segments. Unfortunately, there are few FTR-tip, transient 
response data sets with the above quality. Regardless of the quality of the data, 
however, the least squares method demonstrated itself to be inadequate for even the 
viscous damping only case. As can be seen from figures I4(a)-(f), it has the ten- 
dency to radically overshoot even those portions of the free-tip response curves not 
significantly affected by corner effects. Therefore, the coefficient estimates of 
the logarithmic decay averaging method will be primarily used for the following 
discussion. 
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K contributes only about 51 of the Kg value. K A is approximately equal 
to Kg. Because of difficulty in defining 8 and c, as applied to equation (12), 
the analytical prediction is compared to the experimental Kg estimates. 


There was considerable difference between the error values for each estimate 
of Kg and C fl . Therefore, in comparing the experimental coefficients with the 
predictions from equations (19) and (23) it was advantageous to express the experi- 
mental results in terms of a weighted average that took into account relative error 
estimates. Given the error parameter, Rg^, as defined by equations (8) 
through (10)', a weighting function, W^, can be defined: 


W. 

l 




( 28 ) 


where is the best, i.e., the smallest, envelope-error estimate found in the 

average ensemble. 


Therefore, with this defined function, weighted averaging was performed by the 
following operations: 


N 

£ “A 

K e = ±4 1 (29a) 

where ^ are the spring estimate values. Similarly defined: 


av 


N 

£ 

i=1 


W i C A. 

l 


N 

E 

i=1 


W. 

l 


(29b) 


Having obtained coefficient averages by the above method, an improved compari- 
son can now be made. Figures 15 and 16 are plots of predicted and experimental 
spring and damping constants as a function of rotor rpm. The data includes both 
weighted averages and uncorrected estimates. 

The predicted K ft and C fl values were based on C a estimates derived from 
semispan wind tunnel data, thereby partially taking into account viscous and three- 
dimensional effects. Compressibility effects were not included. However, the 
experimental C a correction provides only approximate predictions of K ft and C A 
for the free tip. The swept-tapered tip C% a used in the C ft and K^ calculations 
was 3.67 rad -1 . The nth harmonic used in the predicted aerodynamic spring and 
damping coefficient calculations was n = 6/rev. This value was chosen because it 
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agreed with the experimental transient response frequency. The selection of n was 
not critical, as and were not sensitive to the harmonic order chosen. 

Establishing the validity of the linear and parabolic trends of the and 
predictions (with respect to V^) are impossible given the limited amount, and 
scatter, of the experimental data. The experimental aerodynamic damping 
coefficients do demonstrate some agreement with the predictions. Almost all data 
points fall within the uncertainty limits of the analytical predictions. These 
uncertainty limits are calculated by the analysis in Appendix C (for limits see 

Appendix D) and the parameter accuracies in table 4. The experimental values, 

though, significantly disagree with the spring predictions. However, a number of 
points do fall within, or on, the uncertainty limits of the predicted aerodynamic 
spring. It is apparent, therefore, that final verification of the free-tip response 
analysis requires an improved data base in terms of both quality and quantity. 

Finally, figures 17(a)-(c) compares the transient response behavior of the free 
tip given the analytical- and experimental-weighted-average dynamic coefficients for 
rotor speeds of V-j, = 1 6 1 , 198, 213 m/sec. The experimental tip time histories were 
not used for the final response comparisons. The nonlinear distortions of the 
experimental time histories make them inappropriate for comparison with the analyti- 
cal predictions. The curves shown in figures 17(a)-(c) were generated by using one 
form of the free-tip response equation's analytical solution 

~^d a)t 

0 = e cos(u).t) 

o d 

where 0 O is the tip's deflection after encountering the air jet. 


SUMMARY 


A motel FTR whirl test was conducted in a small, rectangular test chamber. The 
16-ft diam model FTR underwent collective and rpm sweeps, with and without air- jet 
excitation of the free tips. Both rotor balance loads and tip, pitch-angle time 
histories were acquired. Unlike a previous wind-tunnel test of the FTR, significant 
tip oscillations were observed as a consequence of air- jet excitation. The free-tip 
transient response, with digital filtering, was found to be an approximate single 
degree of freedom, constant coefficient, viscously underdamped oscillatory system. 
However, some nonlinear behavior was observed distorting the transient response. 

This distortion was probably due to test-chamber corner aerodynamic-interference 
effects. 

The results of this test are summarized below: 

1. The FTR free tips responded to the air jet with an approximately 6/rev 
natural frequency. 
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2. Sufficient tip-response data was acquired to estimate free-tip aerodynamic 
spring and damping coefficients. Digital filtering and selective curve fitting were 
necessary to calculate the coefficients. The transient response curves only roughly 
approximate the theoretical exponentially decayed sinusoid. This was possibly due 
to the poor quality of the aerodynamic environment of the test chamber. 

3. Using Theodorsen unsteady, potential-flow analysis analytical expressions 
for the tip aerodynamic spring and damping were developed, K ft and C ft . 

4. Comparisons were made between the experimental and analytical C ft and K fl 
estimates. The aerodynamic spring values showed large disagreements. The damping 
values demonstrated fair agreement. 

5. An analytical expression was developed that related tip geometry and aero- 
dynamics to tip oscillation amplitude, subject to a forcing function having a fre- 
quency that is dependent upon a first order polynomial of J2. It was also analyti- 
cally shown that in hover FTR tip pitching resonance is independent of tip speed. 

6. Limited whirl-test rotor performance data was acquired for low Cp/o. 
Baseline rigid-tip data was not acquired. However, rotor performance data for two 
controller control moments was obtained. No significant performance differences 
were noted for the two controller output settings. 

7. Another free-tip whirl test should be conducted to provide a higher-quality 
data set. Only with such a data set will the analytical K fl and models be 
completely verified. 
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APPENDIX A 


DIFFERENTIAL EQUATION LEAST SQUARES REGRESSION METHOD 


In order to estimate the spring and damping coefficients a curve fit of the 
experimental data could be performed using the tip-response differential equation. 
A least-squares error function was defined by 


E 



+ K E 9 i 


+ sgn(9 i )T f 



(A1) 


where 



10 . 

1 


is the experimental tip inertial moment. The remaining error-function terms com- 
prise the predicted inertial moment. The functional form of E is 

E = E(C A ,K £ ,T f ) 

The error function is subject to the minimization constraint 


where x^ represents the coefficients to be determined by regression, i.e., C A , Kg, 
and Tg. The above constraint results in the following set of equations 


3E 

31 


3E 

3C a 




N 

2 ^ F.e. = 0 
i=1 
N 

2 £ F 1 5 1 = 0 
i=1 
N 

2 ^ F.e. = 0 
i=1 
N 

2^ F.sgn(0.) = 0 
i= 1 


( A2a) 


(A2b) 


( A2c) 


( A2d) 
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where for simplicity is defined to be 

F i = °A®i + K E 9 i + sgn(0 i )T f - M. (A2e) 

1 

The previous set of equations were used to define matrix equations for given 
damping assumptions: viscous damping only, coulomb damping only, and combined 

coulomb viscous damping. These assumptions could be implemented by the selective 
exclusion of constraints and parameters. Experimental data sets were used to calcu- 
late the matrix elements which then were solved by a Gauss-Jordon routine. The 
solution matrix contained the free-tip dynamic coefficients. 

The experimental tip acceleration (e^) was approximated numerically from the 
experimental tip, pitch-angle time-history. Numerical derivatives are inaccurate 
when applied to scattered data, even when incorporating higher-order terms. Yet, it 
was hoped that the data scatter would be minimal and that the constraint summations 
comprising the matrix elements would have a smoothing effect that would compensate 
for the derivative inaccuracy. 

An additional problem with this least squares method is that it is extremely 
sensitive to forced-tip response behavior, such as that induced by test-chamber 
corner effects, and thereby providing inaccurate results. To compensate for the 
test-chamber aerodynamic influence on the response curve, all data-reduction methods 
had the capacity for excluding the forced-tip portion of the transient response. 
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APPENDIX B 


LINEAR DECAY AVERAGING METHOD 


The parameters Kg and T^> are evaluated by this method, which assumes coulomb 
damping only. 

If the tip oscillation damping was completely due to friction, the peak-to-peak 
amplitude decay would be linear according to the relationship 


where 6^ and 6^ + i are adjacent tip oscillation peak amplitudes for each cycle. 

Another consequence of assuming pure coulomb damping is that the tip-response 
frequency is unaffected by the amount of damping. Therefore, 

2ir 

O) = OJ = — 

d x 

and the tip response frequency would be the natural frequency. 

On a local basis, then, the parameters T^ and Kg can be estimated by the 
following equations. 


2tt 



( B1 ) 
( B2) 



9 i + 1 



(B3) 


Finally, averaging is used to calculate a representative set of tip response 
coefficients . 

More extensive discussion of coulomb-damped transient response can be found in 
reference 1 1 . 
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APPENDIX C 


AERODYNAMIC-DAMPING-COEFFICIENT UNCERTAINTY ANALYSIS 


The aerodynamic damping coefficient was found to be 


C A * 5 » 1RC t 
n 


c He - c Vi - c V V T + - c2 (i - c V RV T + 1 ° 3 (i - c Vr 

oyt payyi pay n T n o\4 pay n T 4 o^4 pa J r 


Ca has the functional form 
H n 


C A = C A ( p » AR » C j, » V T» c o’ c pa’ R,F n ,C V 
n n a 


By letting 
C 


'1 


= Ml - c V# - C \f V + - c 2 (l - c \g RV t + l c 3 (j - c VJ 
o\4 paM4 pay n T n ol4 pa J n T 4 oy4 pay Tj 

» t » 


C 2 = g P ARC^ V T 
a 


The following partial derivatives can be defined 


3G A 1 

IT - 5 4RC i Vi 

a 


3G A 1 

iZR = 2 pC l V T C 1 


(Cla) 


(Clb) 


3G A 1 

icT = ^ p arv t c i 


(Clc) 


3G A 1 

Iv! = 2 p ARC * c i 

T a 


^A 

3C = °2 
o 


3C a 

— — - c 

3c " 


pa 


\sc 2 ljr - c Vj -c )f + -c(j-c'\gR + J c 2 (f - c Y| 
o\4 pay^4 pa J n n o^4 pa ) n 4 oyJ payj 

|*c 3 F (2c - 1) - - G Rc 2 - l c 3 l 

o n pa n n o 4 ol 


(Old) 


(Cle) 


(Clf) 
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(Clg) 


( C 1 h ) 


(Cl i) 


The uncertainty in F n and G n is governed by the uncertainty in k, the 
reduced frequency. Noting 


k 


c nfi 
o 



c n 
o 


" 2R 


where functionally 


k = k(c Q ,R) 

Its partial derivatives are therefore 

3k _ n_ 

3c = 2R 
o 


(C2a) 


IK 

3R = 


c n 
o 


2R 


(C2b) 


A computer algorithm was developed to perform the uncertainty calcula- 

tions. The uncertainty analysis is based upon a method presented in reference 12. 
This method requires the definition of the previous derivatives to define u>£^, the 
uncertainty of the predicted C^. One simplifying approximation used in the C A 
uncertainty calculation is the uncertainty estimate for Y i ( k ) , Y Q (k): only the 

accuracy of the reduced frequency is considered. The accuracy of the Bessel func- 
tions J Q (k) and J i ( k ) influencing the accuracy of Y Q (k) and Y^(k) is not incorpor- 
ated into the uncertainty analysis. 

The uncertainty of F n and G n can now be estimated once the Bessel function 
uncertainties have been estimated from the k accuracy estimate. 


G 


n 


j i tY i - v - v j i * v 

or, - y 2 * <j, * y 2 
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where 


Letting 


G = G (J ,Y , Y, ) 
n n o’ o’ 1’ 1 


G1 = J 1< Y 1 - Jo) - Y 1 (J 1 + Y o } 
G2 = (Y, - J _) 2 + (J. + Y ) 2 


The G n derivatives are 


3G n 

- J 1 

2G1(Y 1 - J Q ) 

3 J 

0 

' G2 

+ 2 
G2 

3G 

n 

" Y 1 

2G 1 ( J 1 + Y q ) 

3Y 

0 

* G2 

2 

G2^ 

3G 

n 

-J 

o 

2G1(J 1 + Y q ) 

3J 1 

= G2 

2 

G2 

3G 

n 

-Y 

o 

2G 1 ( Y 1 - J q ) 

3Y 1 

= G2 

2 

C2^ 


(C3a) 

(C3b) 

(C3c) 

(C3d) 


Similarly for F n 


where 


Letting 


F J 1 (J 1 * , o> - *1»1 - J c> 

" ’ - J o » 2 * < J , * V 2 

F n * F n (J o’ Y o' J 1-V 

p 1 • J,(J, * 1 g ) ♦ Y,U, - J c ) 

F2 = <Y, - J 0 ) 2 * (J, * y/ 
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partial derivatives are 


The F n 


3F 

n 

- Y 1 

2F1(Y 1 - J q ) 

3 J 

0 

" F2 

+ 2 
F2 

3F 

n 

J 1 

2F1 ( J. + Y ) 
1 0 

3Y o 

“ F2 

F2 2 

3F 

n 

(2J 1 

. Y 0 ) 2F1(J, . Y^) 

3 J 1 


F2 - F2 2 

3F 

n 

(2Y 1 

- J o> 2F1 <*1 * V 

3Y 1 


F2 F2 2 


(C4a) 


(C4b) 


(C4c) 


(C4d) 


The final step in defining the uncertainty of the predicted C A 
can be accomplished by the following general expression (ref. 12). 


results, 



1/2 


where x^^ represents all the independent parameters, is the estimated accuracy 

of each independent parameter, and wq is the uncertainty of the dependent (calcu- 
lated) parameter, C. 


Therefore, for 


C A the previous expression becomes 



(C5) 


The predicted aerodynamic damping uncertainty can be calculated by first using 
the parameter values and accuracies in table 4 to calculate the partial derivatives, 
equations (C1)-(C4). Equation (C5) is then used to calculate u)c ft . For the swept 
and tapered free tip used in this test, u>c A was typically 411& of the predicted 
C A values. 
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APPENDIX D 


AERODYNAMIC SPRING UNCERTAINTY ANALYSIS 


The aerodynamic spring equation was found to be 


k a = 5 » 4RV ?°o 
n 


(s - °pa)* C A a [ F n - C c(? - c pa) ¥] 


Therefore, in functional form 

K 

By letting 


A = K A ( p ’ AR ’ V T ,c o ,c pa’ C Jl ,F n’ G n ,R ) 
n n\ K a / 
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1 o^ pa kj n oyi pay R 


K 2 = \ p ARV 2 

The partial derivatives are then 
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= K 2 c 


o(°pa * ?)(! - C pa) C l a („2 ) 


ar = k 2= 

n 


c(°pa - 4 ) C l a 


3G* = K 2°c(° P a ' 4)^(4 - 0 pa)(s) 


(Dlh) 


(DU) 


Finally, to estimate 3K A /3F n and 3K A /3G n the contribution of k, the reduced 
frequency, must be incorporated into the uncertainty analysis (refer to Appendix C). 

The uncertainty analysis method of reference 12 is used to determine the uncer- 
tainty of the predicted results, i.e., ui[( A . Using table 4 parameter values and 

accuracies the K A partial derivatives can be calculated and substituted into an 
uncertainty estimate equation similar to equation C5. Typical values are 39# 

of the predicted values of the free tips that are whirl-tested. 
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TABLE 1.- ROTOR AND FREE-TIP GEOMETRY 


Model rotor (dimensions include free tips) 

Radius 2.56 m 

Number of blades 4 

Chord 0.171 m (constant to 0.9R) 

Solidity 0.085 

Twist (center of rotation to tip) -9.45° (for A9 = 0) 

Airfoils V23010-1.58 (constant) 

Cutout 0. 1825R 

Flap-hinge offset 0.031R 

Weight moment about flap hinge 46.8 N- 

Inertia about flap hinge 6.17 kg - m 


Free-tip geometry ( 10% tip) 

Extends over outboard 10% of rotor radius 

0.30 taper ratio over outboard 5% of blade radius 

35° of 1/4 chord sweep over outboard 5% of blade radius 

V23010 - 1.58 with 0° tab angle 

Tip twist, -0.94° 

Pitch axis, 13% of chord 

Center of gravity located at 14% chord 
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TABLE 2.- HARMONIC COEFFICIENTS FOR TEST DATA POINTS 


Whirl-test rotor rpm 

796 

740 

600 

Test data point 

1 

2 

3 

4 

5 

6 

Coefficient 







A 1 

0.212 

0.334 

0.485 

0.212 

0.371 

0.408 

a 2 

-.019 

-.051 

-.087 

-.018 

-.063 

-.121 

a 3 

.062 

.097 

.079 

.054 

.180 

. 181 

a 4 

-.137 

-.176 

-.188 

-.085 

-.127 

-.284 

A 5 

.069 

.175 

.193 

.053 

.141 

.291 

a 6 

-.015 

-.077 

-.112 

-.022 

-.112 

-.247 

a 7 

.025 

.066 

.096 

.015 

.069 

.174 

a 8 

.000 

-.027 

-.099 

-.030 

-.055 

-.100 

a 9 

-.016 


-.002 

.016 

.033 

.047 

A 10 

-.008 

9 

CTn 

o 

o 

1 

-.002 

-.011 

.014 

B 1 

.057 

.161 

.268 

.061 

.086 

.245 

b 2 

-.067 

-.095 

-.056 

.001 

-.064 

-.184 

b 3 

.055 

.045 

.113 

.047 

.061 

.111 

b 4 

-.029 

-.065 

-.099 

-.014 

-.057 

-.142 

b 5 

-.056 

-.018 

.044 

-.009 

.005 

.063 

b 6 

.048 

.069 

.064 

.018 

.076 

.125 

b 7 

-.033 

-.070 

-.051 

-.015 

-.046 

-.137 

b 8 

.039 

.093 

.064 

.018 

.059 

. 132 

B 9 

-.023 

-.023 

-.117 

-.007 

-.050 

-.133 

B 10 

.018 

.042 

.053 

.016 

.071 

.158 


The averaged free-tip response can be represented by a Fourier 
series 

10 

6 = ^ [A n sin(nt) + B n cos(nt)] 
n=1 

where Ap A^q and B^, B^q can be obtained 

from the above table for six test points and three different 
rotor rpm's. 
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TABLE 3.- DATA-REDUCJION ALGORITHMS MATRIX 


Data-reduction 

method 

Calculates 

free-tip 

spring 

Calculates 

free-tip 

viscous- 

damping 

Calculates 

free-tip 

coloumb- 

damping 

Differential-equation least-squares fit 

Yes 

Yes 

Yes 

(complete coefficient set) 




Differential-equation least-squares fit 

Yes 

Yes 

No 

(partial coefficient set) 




Logarithmic-decrement averaging 

Yes 

Yes 

No 

Linear-decay averaging 

Yes 

No 

Yes 


TABLE 4.- FREE-TIP AERODYNAMIC SPRING 


AND 

DAMPING PARAMETERS AND 

ACCURACY 

P 

1 .227x10“-’ kg/m 3 

±1 .227x10 -6 

AR 

0.2097 m 

±0.0210 

c aa 

0.064 deg -1 

±0.0064 

C o 

0. 1710 m 

±0.0171 

C pa 

0.13/5 chord 

±0.013 

V t 

188.4578 mps 

±2.5603 

R 

2.4384 m 

±0.0610 
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d? POOR QUALITY 


PITCH 

AXIS 


CONTROL 

MOMENT 


SELF ADJUSTING 
BLADE PITCH 


PITCH AXIS AND 
CENTER OF GRAVITY 





Figure 4.- Controller and rotor tip assembly. 



Figure 5.- Free-tip jet response time history. 
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kPa 

137.90 


•-A-* 413.69 
— O- 275.79 



LONGITUDINAL MEASURING STATION LOCATION 
OF THE RECTANGULAR JET, m 

(a) Jet-velocity profile in the rotor tip path direction (vertical distance from 

the jet = 0.46 m). 



LONGITUDINAL MEASURING STATION LOCATION 
OF THE RECTANGULAR JET, m 


(b) Jet-velocity profile in the rotor tip path direction (vertical distance from 

the jet = 0.76 m) . 

Figure 6.- Measured air-jet velocity distributions. 
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LATERAL MEASURING STATION LOCATION OF 
THE RECTANGULAR JET, m 


(c) Jet-velocity profile in the rotor spanwise direction (vertical distance from 

the jet = 1 .83 m) . 

Figure 6.- Concluded. 
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180 220 260 300 340 380 

AZIMUTH, deg 


Figure 7.- Experimental free-tip Jet-response time history, digitally filtered. 
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MOMENT, N-m 



0 20 40 60 80 100 

TWIST ANGLE, deg 


(a) TSC moment vs. twist characteristics. 


5r 

--0-- CALIBRATION POINTS 



0 1500 3000 4500 6000 7500 


RADIAL LOAD, N 

(b) TSC moment vs. tensile load characteristics. 

Figure 10.- Free-tip torsion-strap controller pull test results. 







(c) M) 


Figure 11.- The envelope-area error estimate, (a) Transient response with envelope 
outlined, (b) The envelope area, (c) Curve-fitted response and experimental 
data, (d) The resulting envelope area comparison. 
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TOP VIEW 



CONTROLLER OUTPUT 
SHAFT BEARING 
SURFACES 



CONTROLLER TENSILE 
FORCE = CF cosj3 cos?* 


SIDE VIEW 


CONTROLLER/FREE 
TIP C.G. 



Figure 12.- A freebody diagram of the force components contributing to the 

controller straps' tensile load. 



TOP VIEW 



CF cos P sin ?* 


SIDE VIEW 


CONTROLLER/FREE 
TIP C.G. 


OUTPUT SHAFT / 

BEARING SURFACES / CFcos/3 



Figure 13.- A freebody diagram of the force components contributing to the 

controller bearings' contact force. 
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Figure 14 



EXPERIMENTAL (796 rpm, TP = 4) 



(b) Test point 2 averaged transient response, 796 rpm. 

The filtered free-tip transient responses used for the spring and 
damping estimates. (Curve-fit results also shown.) 
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EXPERIMENTAL (740 rpm, TP= 1) 
LOGARITHMIC DECREMENT 
LEAST SQUARES 



AZIMUTH, deg 


(c) Test point 3 averaged transient response, 740 

020 1" EXPERIMENTAL (740 rpm, TP = 2) 

.015 - LOGARITHMIC DECREMENT 

LEAST SQUARES 


AZIMUTH, deg 


(d) Test point 4 averaged transient response, 740 


Figure 14.- Continued. 



6, rad 0, rad 


EXPERIMENTAL (600 rpm, TP = 1) 



Test point 5 averaged transient response, 600 rpm. 


EXPERIMENTAL (600 rpm, TP = 2) 

LOGARITHMIC DECREMENT 

LEAST SQUARES 



(f) Test point 6 averaged transient response, 600 rpm. 
Figure 14.- Concluded. 
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O EXPERIMENTAL 
□ WEIGHTED EXP. AVERAGES 

ANALYTICAL PREDICTIONS 

150r PRED. UNCERTAINTY LIMIT 

O O 


s 

/ 

/ 



400 600 800 1000 

ROTOR rpm 

Figure 15.- Free-tip aerodynamic-spring analytical and experimental comparison. 


O EXPERIMENTAL 

□ WEIGHTED EXP. AVERAGES 

ANALYTICAL PREDICTION 

.15- PRED. UNCERTAINTY LIMIT 



400 600 800 1000 
ROTOR, rpm 

Figure 16.- Free-tip aerodynamic-damping analytical and experimental comparison. 
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TIME, sec 



TIME, sec 

Figure 17.- Analytical and empirically generated transient responses. 

(a) Transient-response comparison, 796 rpm. (b) Transient-response comparison, 
740 rpm. (c) Transient-response comparison, 600 rpm. 
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